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ABSTRACT 

Numerical relativity simulations predict that coalescence of supermassive 
black hole (SMBH) binaries not only leads to a spin flip but also to a recoil- 
ing of the merger remnant SMBHs. In the literature, X-shaped radio sources 
are popularly suggested to be candidates for SMBH mergers with spin flip of 
jet-ejecting SMBHs. Here we investigate the spectral and spatial observational 
signatures of the recoiling SMBHs in radio sources undergoing black hole spin 
flip. Our results show that SMBHs in most spin-flip radio sources have mass 
ratio q > 0.3 with a minimum possible value gmm — 0.05. For major mergers, the 
remnant SMBHs can get a kick velocity as high as 2100 km s~^ in the direction 
within an angle < 40° relative to the spin axes of remnant SMBHs, implying that 
recoiling quasars are biased to be with high Doppler-shifted broad emission lines 
while recoiling radio galaxies are biased to large apparent spatial off-center dis- 
placements. We also calculate the distribution functions of line-of-sight velocity 
and apparent spatial off-center for spin-flip radio sources with different apparent 
jet reorientation angles. Our results show that the larger the apparent jet reori- 
entation angle is, the larger the Doppler-shifting recoiling velocity and apparent 
spatial off-center displacement will be. We investigate the effects of recoiling 
velocity on the dust torus in spin-flip radio sources and suggest that recoiling 
of SMBHs would lead to "dust poor" AGNs. Finally, we collect a sample of 19 
X-shaped radio objects and for each object give the probability of detecting the 
predicted signatures of recoiling SMBH. 

Subject headings: black hole physics — galaxies: active — galaxies: evolution — 
galaxies: jets — gravitational waves — quasars: general 



Introduction 



In the hierarchical galaxy formation mod el of the ACDM cosmology, massive galax- 
ies are produced by successive galaxy mergers (IKauffmann &: Haehneltll2000l : ISpringel et al, 
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20051 ). During the merging of galaxies, the supermassive black holes (SMBHs) at centers 
quickly form hard supermassive black hole binaries (SMBHBs) at a separation of j)c-scale, 
and then may stall for a timescale longer than Hubble time, if the gravitational potential 



(Beeelman. et al. 


1980; 


Yu 


2002) 



and SMBHBs become merged within Hubble time because of triaxial or rotating structures 



massive perturbers, or massive gas disks at galactic centers ( iMerritt fc Milosavljevic 112005 



Colpi fc Dottill201ll . and references therein). The strong gravitational wave (GW) radiations 
from the coalescence of SMBHBs are the main targets of the GW detector Laser Interferom- 
eter Space Antenna (LISA) and the GW detecting program Pulsar Timing Array (PTA). 

Spatially resolved dual or binary active galactic nuclei (AGNs) at fcpc- scale separation 



2006; 


Green et al. 20 IC 


1 ) , and recently in a large sample 


of AGNs showing d 


ouble-peak narrow 


emission lines 


Zhou et al. 2004: Xu & Komossa 2009 


; IComerford et al. 


2009; IwauET et al. 


2009; 


Liu et a 




2010a 


ISmith et al. 201ol: Liu et al. 2010b; 


Fu et al. 2010 


Shen et al. 2011; 


Rosario et al. 


201ll). Meanwhile, more compact or hard SMBHBs are extremely difficult to 



observe with present telescopes (see [Rodriguez et al.1 l2006l for the only pc-scale SMBHB 
candidate so far discovered by VLBA), and for sub-pc SMBHBs only indirect observa- 
tional evidences are su ggested in a few sub-classes of AGNs with peculiar ra dio morphologies 



(iBege. 



1988 



man, et al 



1980 



Valtaoja et al 



2000 



emission-line features ( 



Liu et al 



1995 



1997 



Stockton fc Farnham 



Liu 1120041: iLiu fc Chen II2007I). periodic variabilitie s (jSillanpaa et al. 



2006 



1991 



Liu fc Wu Il2002[). or special broad - 



Gaskell I Il996l : iBoroson fc Laueii 120091 ). 



Although it is still unclear how SMBHBs can pass through the final pc-scale and become 
merged, some observations of AGNs and QSOs are suggested to be due to the coalescence of 



(Liu et al. 


2OO3I) 


2002 


: Izier 


2005) 



20031) . and the rapid jet reorientation in X-shaped radio sources (IMerritt &: EkersI 



20051 ) . For SMBHBs i n quiescent galaxies some observationa l signatures have 



been predicted ( IChen et all I2OO8L I2009L 1201 ll : iKomossa fc MerrittI l2008bl : iLiu et all |2009|) 



but need to be confirmed observationally. Among all the observational evidences for the 
coalescence of SMBHBs in AGNs, X-shaped radio source is one of the most popular. X- 
shaped radio sources are a subclass of double-lobed radio sources, consisting of a pair of 
l ong and diffusive inactive "wings" oriented at a la rge angle relative to the active lobes 
( iLeahy fc Williamslll984t iDennett-Thorpe et al.ll2002l ). The "spin fiip model" suggests that 
the wings are the fossil plasma of old lobes and jets are misaligned with the wings because 
of swift jet re orientation due to the spin-fiip of jet-ejecting SMBHs during the coalescence 
of SMBHBs ( iMerritt fc EkersI |2002| ). Meanwhile, many other formation scenarios, which 
do not require binary black hole coalescence, are also proposed in the literature, including 
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merger of two jet-emitting AGNs 
and intra- /inter-galactic medium I 


(Lai & Rao 


2007) 


, hydrodynamic interaction between iet 


Wirth et al 


. 1982 


: Leahy & Williams 


1984: ^ 


A'^orrall et al. 


1995 


CaDetti et al. 2002: Zier 2005: SariDalli & Subrahmanvan 2009: 


Hodees-Kluck et al. 


2010 


, and precession of jet axis (] 


5kers et al. 


19781; 


Beeelman. et al. 1980: Liu 


2004). 



Numerical simulations show that the coalescence of two black holes not only leads 
to a spin-flip black hole but also to a kick of the post-merger black holes due to the 
anisotropic emissions of gravitational waves. Therefore, it is expected that the coales- 
cence of SMBHBs in spin-flip radio sources should lead to the recoiling of post-merger 
SMBHs surrounded by accretion disk. The recoiling velocity depends on the angular mo- 
mentum and mass ratio of th e pre-merger SMBHs , and can be as large as 175 km s~^ 
for Schwarzschild black holes (IGonzalez et al.l 120071 ). 3800 km s~^ for extremely spinning 



( Pretorius 


2005: 


2007a 


b), or even 



Herrmann et al.l 120071 : iBaker et al.l |2006|. 



b|), or even lO^kms ^ in hyperbolic encounters (IHealy et al.l l2009l ). The gravita- 



2008 



Campanelli et al.l 12006 



tional r ecoiling SMBHs in AGNs may lead to s i gnatures such as kinematic and spatial 



offsets (|Bonning et a 



Robinson et a 



Magain et al 



2010 



200 



i 



2007 



Komossa et al. 



Tsalmantza et al. 



Loebll2007l: 



2011 



2008 



Shields et al. 



Eracleous et al.l 



Batcheldor et al 



20101: 



2009b 



2011 



Civano et al 



Madau fc Quataert 



Jonker et al 



201 



3 



2010 



2004 



Blecha fc Loeb 



20081 : iGuedes et al.ll2009l : ICivano et all 12010 



Guedes et a 



with some other electromagnetic signatures (see 



Komossa 



20 111 : iBlecha et al.ll201l[ ). along 



2010l . for a review) 



In the following, we refer to the radio sources with spin-flip signatures, such as those 
in X-shaped radio sources, as "spin-flip radio sources" . In this paper, we adopt the spin-flip 
model and investigate the spectral and spatial signatures of the recoiling SMBHs in such 
objects, using the empirical formula obtained in numerical relativity to calculate the spin 
angular momentum and recoiling velocity of the post-merger black hole. We neglect the slow 
precession of the spin axes of SMBHs during the in-spiraling of SMBHBs, We assume that 
the old and active jets, respectively, form along the spin axes of the primary and post-merger 
SMBHs and any jet-ejecting SMBH has spin parameter a = cS/GM"^ > 0.9, where S is the 
black hole spin angular momentum, M is the black hole mass, G is the gravitational constant, 
and c is the speed of light. For the secondary SMBH, no constraint on the spin parameter 
a2 can be given. We assume that 02 can be any value between and 1 of flat probability 
distribution. Our results show that the SMBHBs producing detectable spin-flip angles should 
form in major mergers with mass ratio > 0.2—0.3 and that the recoiling SMBHs preferentially 
move in a direction with an angle < 40° relative to jet orientation. We calculate the Doppler 
kinematic shift of broad emission lines, the AGN off-center displacement in host galaxy, and 
the effects of recoiling SMBHs on dust torus in spin-flip radio galaxies and quasars as a 
function of apparent jet reorientation angle and jet size. We finally collect a sample of 19 X- 
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shaped radio objects and calculate the detection probabilities for the proposed observational 
signatures. 

The paper is organized as follows. In Section [21 we introduce the empirical formula 
for the spin angular momenta and the recoiling velocity of post-merger SMBHs obtained in 
numerical relativity. Section [3] describes our numerical simulations and general results for 
spin-flip radio sources. We collect a sample of 19 X-shaped radio objects and calculated the 
observational signatures of recoiling SMBHs in Section HJ Our discussions and conclusions are 
given in Section |51 Throughout the paper, we assume a ACDM cosmology with parameters 
= 73 km s-i Mpc~\ Qa = 0.73, and = 0.27. 



2. Empirical formula for spin angular momenta and recoiling velocities of 

post-merger SMBHs 

It is shown that the final spin and the recoiling velocity of a post-merger SMBH can 
be calculated in high precision with the emp i rical f ormula fitting the simulation results of 
numerical relativity. Following iRezzoUa et al.l ( l2008l ). the final spin angular momentum Sfin 
of a post-merger SMBH can be approximated as the sum of the two initial spin angular 
momenta Si and S2 of the primary and secondary SMBHs, respectively, and a third vector 



'fin 



S1 + S2 



where i is the difference between the orbital angular momentum L at large SMBHB separa- 
tion and the angular momentum Jrad radiated away during merger 



L- J 



rad- 



(2) 



Equation ([T]) can be written as 

3-fin 



(ai + asg^ + iq) 



(3) 



(iRezzoUa et al.l 120081 ) . where mass ratio q = M2/M1 < 1, spin parameter vector agn 



Sfin/M^, £ = £/{MiM2), the primary SMBH spin parameter vector ai = Si/Mi^, and the 
secondary SMBH spin parameter vector a2 = S2/M2^. M, Mi, and M2 are, respectively, 
the masses of the final SMBH, the primary SMBH, and the secondary SMBH. Here and in 
section (j2]), we use the natural units G = c = 1. In equation ([3]), the direction of the angular 
momentum vector i is nearly parallel to the orbital angular momentum and its norm is 



1^1 



S4 



+ |a2|^g'' + 2|aip|a2|^g^ cos a) + 
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S5V + to + 2 

+ t2r] + tsT]^ 



^ (|ai| cos/3 + |a2|g^ C0S7) + 



(4) 



where S4, S5, to, h, and t-s are the empirical fitting coefficients given by lRezzolla et al.l ( 120081 ) . 
rj is the symmetric mass ratio rj = M1M2/ {Mi + M2Y = 9/(1 + qYi and the three projected 
(cosine) angles a, (3 and 7 are defined with the inner products 



cos a = ai ■ 3.2 , cos (3 = ki ■ 



cos 7 = a2 



(5) 



In the definition of the orbital angular momentum L is the one at large (infinity) separation. 
In the calculation, we neglect the slow precession of spin axis around total angular momentum 
and use all the empirical formula for a finite separation ^ M to calculate the spin-fiip 
angle. We expected that this approximation can give good enough results to describe the 
jet reorientations in X-shaped radio sources. 

Coalescence of SMBHBs not only leads to the spin-fiip but also to a fast recoiling of the 
post-merger SMBH. The recoiling velocity of post-merger SMBH can be empirically given 

by 



V 



Vm ei + V^s (cos ^ ei + sin { 62) 
Ar/Vl -4r/(l + 5r/), 



H- 



7] 



.11 



(6) 



(7) 



( ICampanelli et al.ll2007a] : iGonzalez et al.l 120071 ). where ei and 63 are, respectively, the unit 
vectors in the directions of the binary separation from the primary to the secondary black 
hole and the orbital rotating axis just before merger, and the unit vector 62 = ei x 63. The 
index || and _L refer to, respectively, the projections of the vectors parallel and perpendicular 
to the orbital axis. The index m indicates the recoil velocity of unequal mass contribution, 
and s indicates the contribution due to spin. The fitting parameter ^ is the angle between 
the recoiling velocities in the orbital plane due to the mass a symmetry and spin a ngula r 
momenta and we use ^ = 145° in the following calcul ations as inlLousto fc Zlochowerl (120081 ). 
The fitting coefficients A, B, and H are taken from ICampanelli et al.l (l2007al). For the kic k 
velocity parallel to orbital axis, we use the empiric al formula given bvlyan Meter et al.l (120101 ). 
These formula are equivalent to one suggested by lLousto et al.l (l2010al ) if angular parameters 
are suitably interpreted. In order to use these forrn u la for statistical calculation, we combined 
the equations (7), (8) an d (9) in Ivan Meter et al.l (120101 ) and adopted the interpretation of 



angular parameters from lLousto et al. 



(l2010al ). Then we obtained the following formula 



— I cos(6a — ©o) 
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Ksjq - IW 



+ g 02 I cos(6s' - 61^ 



(8) 



where K2, K^, and Kg are fitting parameters taken from Ivan Meter et al.l (120101 ). In equa- 
tion ([8]), af- is tlie magnitude of the spin parameter vectors ai of the ith black hole pro- 
jected into the orbital plane, 6a (©5) is the angle between the in-plane component of 
A = S2/M2 — Si/Mi (S = Si -|- S2) and the unit separation vector from the primary toward 
the secondary at merger and Gq (Gi) is the corresponding initial value of 6a (©s) at the 
arbitrary separation of the binary with quasi- circular orbit. 



3. Numerical Simulations and Results 

3.1. Monte Carlo simulations 

With the empirical formula given in Section [21 we now investigate the distributions of 
the SMBH mass ratio, spin-flip angle, recoiling velocity in spin-flip radio sources under the 
spin-flip model. In our calculations, the physical quantities are specified in the Cartesian 
coordinate system defined by ei, e2, and 63 and so that the positive direction of the 2;- axis 
align with the orientation of 63 and the orbital plane of the pre-merger SMBHB is in the x-y 
plane. To calculate the final spin vector afin and the recoiling velocity V for a post-merger 
SMBH, we have to specify the following physical quantities of the pre-merger SMBHB: q, 
©A, ©5, ©0, ©1, and both the magnitude and orientation of ai and a2. We employ the Monte 
Carlo scheme to generate these quantities according to the statistical properties motivated 
by the physical considerations. The details follows. 

The mass ratio q of the SMBHB has a uniform distribution in the range [0,1], and 
the spin vec t ors ai and a2 are randomly oriented (unless otherwise noted). Following 
Lousto et al.l f l2010a| ) and for statistical purpose, in the simulations we define ©a and ©s 



with respect to ei, and take ©0 = 0, ©1 has a uniform distribution in the range [0,27r]. 
In the case of an spin-flip radio source, although both SMBHs before merger could be the 
candidates for the jet-ejecting black holes, we assume that the low luminosity diffusive lobes 
are produced by the primary SMBH before coalescence and the active jets form along the 
final spin axis of the post-merger SMBH. Because we are interested in radio loud AGNs or 
quasars with powerful relativistic jets both before and after the merger, we do simulations 
only for uniform distributions of Oi within the range 0.9 < Oi < 0.998. For the spin of 
the secondary SMBH, we cannot give any priori constraint on it thus we assume a uniform 
distributions in the range of < 02 < 0.998. Among the simulations, we adopt only those 
corresponding to 0.9 < afin < 0.998 as the final results for spin-flip radio sources. Although 
we set an upper limit 0.998 rather than the extreme spin a = 1 to all the spin parameters by 
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physical consideration, our results are not significantly changed by this small difference in 
upper limit because of the regularity of the results and the rarity of SMBH merger at a = 1. 

In each run of Monte Carlo simulations, we first draw a set of physical quantities ac- 
cording to the distribution functions specified above, and then calculate agn and V, using 
the equations (I3])-(I8]). To derive the distribution functions of the physical quantities, we 
use 10''-10^ sets of the quantities, depending on the request of resolution. To compare the 
numerical results with the real observations, we project the vectors agn and V both along 
and vertical to the line of sight (LOS) vector eg, which is assumed to randomly orient. We 
calculate the angles wi and W2, respectively, between ai and Sg and between afin and eg, 
which may be observed directly or observationally constrained. Because we assume that jets 
form along the spin axes of the primary SMBHs before merger and of the final SMBHs after 
merger, the spin flip during merger leads to the reorientation of jets. We also calculate the 
projected (apparent) jet reorientation angle 6, i.e. the angle between the components of ai 
and afin projected vertically to LOS 

cos(5) = (^iX^s)-(aHnXe)_ 

I ai X Gg 1 1 agn X Gs I 

To be identified as a spin-flip radio source, both the active jets forming along the spin 
axis of the post-merger SMBH and the relic jets forming along the spin axis of the pre-merger 
primary SMBH should be resolvable by radio telescopes. Moreover, the active jets should 
orient at a large angle with respect to the relic jets. In the special case of X-shaped radio 
sources, the length of the relics jets (wings) should be at least 80% o f the length of the activ e 



ones, and the angle between relic and active jets should exceed 15° ( iLeahy fc Parma Ill992l ). 
The requirements of length and symmetry imply that the angle between relic jets and LOS 
cannot be very small, and clear identification of diffusive extended wings from the plasma 
cocoon of active jets also requires a large relative angle between the projected wings and 
active lobes in the plane of the sky. Therefore, the simulations resulting in wi < 20° or 
6 < 15° are not included in the later statistical analysis. Besides, a radio source would be 
observed as a blazar if the relativistic jet is nearly along the line of sight and the relativistic 
beaming effect becomes dominant. Resolving active jets and good measurement of their 
length and orientation therefore require that W2 cannot be very small. For this reason, the 
simulations resulting in W2 < 10° are excluded. 

After excluding the solutions of the equations which do not fulfill the above selection 
criteria, only about 10^-10^ simulations are left and will be used in the derivation of the 
distribution functions of the SMBH mass ratio q, the spin fiip angle A, and recoiling velocity 
V for different apparent jet reorientation angles 6, and in the investigation of the effects of 
recoiling velocity on dust torus. To analyze the standard deviation of a distribution function. 
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we run in total 5 sets of Monte Carlo simulations. 



3.2. Distributions of SMBH mass ratios in Spin-flip radio sources 

Our results of the distribution of SMBH mass ratio q in the simulated spin-flip radio 
sources are shown in Figure [TJ We calculate the apparent jet reorientation angle using 
Equation ([9]) and divide the simulation results into three groups, based on the apparent 
jet reorientation angles for 15° < 5 < 40°, 40° < 5 < 60°, and 65° < 5 < 90°. Given the 
conditions we have, the mass ratio q in spin-flip radio sources cannot be uniquely determined, 
therefore we are able to give only the statistical probability q{dP/dq) of a certain mass ratio, 
as the distribution function in Figure [H 

Our results imply that to form a detectable apparent jet reorientation angle with S > 15°, 
the mass ratio of the S MBHB must be larger than 0.05, which is smaller than the minimum 



mass ratio obtained by iHughes &: BlandfordI (120031 ). This difference of the minimum mass 



ratio is due to the projection effects on the spin flip which can amplify a small jet reorientation 
angle. If we set smaller lower limits on wi and W2, the amplification effect becomes more 
significant and we may get an even smaller minimum mass ratio. Our results also show that 
for a larger apparent jet reorientation angle the needed minimum mass ratio becomes larger. 
To form spin-flip radio structure with apparent jet reorientation angle 40° < 6 < 60°, the 
needed minimum mass ratio is gmin — 0.1, while for 65° < 6 < 90°, q^i^ — 0.15. 

Figure [1] showed that for an spin-flip radio source with apparent jet reorientation angle 
15° < 6 < 40°, the SMBHB most probably has mass ratio q ~ 0.3, and large contribu- 
tion comes from SMBHBs with mass ratio 0.15 ^ q ^ 0.5 although the contribution from 
SMBHBs with q > 0.5 cannot be completely neglected. For spin-flip radio sources with 
larger apparent jet reorientation angles, the contribution of the mergers of SMBHBs with 
larger mass ratio q become more significant while that of SMBHBs with smaller mass ratio 
becomes less important. For spin-flip radio sources with intermediate apparent jet reorien- 
tation angle 40° < 6 < 60°, the mergers of SMBHBs with mass ratio q > 0.25 may make 
significant contributions, although the contribution from SMBHB mergers with q ~ 0.5 is 
at local minimum. For 65° < 6 < 90°, the jet reorientation in spin-flip radio sources is 
likely due to the merger of SMBHBs with mass ratio q > 0.35, the larger the mass ratio, 
the more likely, and the largest contribution is due to the mergers with q > 0.6. We notice 
that for 65° < d < 90° the typical probability q{dP/dq) could be larger than unit although 
the integrated probability is unit. This is due to the reason that the probability per unit 
mass ratio varies very rapidly on a length scale much less than the typical mass ratio or 
P/{dP/dq) < g for g > 0.8. 
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Fig. 1. — Distribution function of SMBHB mass ratio q in spin-flip radio sources for different 
apparent jet reorientation angles 5. Apparent jet reorientation angle 5 is the angle between 
the orientations of the active and remnant (wings) jets, which is also the projected angle 
between the spin axes of the primary and the post-merger SMBHs. The minimum mass 
ratios depend on 5 and are, respectively, gmin — 0.05 for 15° < 5 < 40°, gmin — 0.1 for 
40° < 5 < 65°, and gmin — 0.15 for 65° < 5 < 90°. SMBHBs in most spin-flip radio sources 
have mass ratio q > 0.3 and thus are formed by major mergers. Error bar gives Icr standard 
error. 
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3.3. Spin flip angles and recoiling velocities of post-merger SMBHs 

Our results suggest that the merger of SMBHB with a given mass ratio can lead to 
apparent jet reorientation angle 5 in a broad range, which is partially due to the projection 
effect. In this section, we show the distribution function of the intrinsic spin flip angle A 
We calculate A using the relation cos(A) = ai ■ afin/|ai||afin| and give the results for typical 
mass ratios q = 0.3, q = 0.7, and g = 1 in the left panel of Figure [2l In the same panel, we 
also give the distribution of A in minor merger g = 0.1, which is also the lower mass ratio 
limit to form spin-flip radio sources with moderate apparent jet reorientation angles. 

We can see that in spin-flip radio sources the typical spin flip angle corresponding to the 
maximum probability increases with mass ratio q from A = 12° for minor merger q = 0.1, 
A = 19° for major merger q = 0.3, and A = 25° for major merger q = 0.7 to A = 30° 
for equal- mass merger q = 1. Our results also show that the spin flip angle A has a broad 
distribution and the full width at half maximum FWHM(A) also increases with SMBHB 
mass ratio from FWHM(A) ~ 7° for minor merger g = 0.1 through FWHM(A) ^ 13° for 
major merger q = 0.3 to FWHM(A) ~ 28° for equal-mass merger q = 1. For equal-mass 
mergers (g = 1), we flnd A < 55°. The trial solutions with A > 55° are mostly rejected 
because that signiflcant spin down of the post-merger SMBH would occur. 

The maximum spin flip angle for minor merger g = 0.1 is about 15°, which is consis- 
tent with the results given in Figure [1] in which no spin-flip radio source with apparent jet 
reorientation angle 6 > 40° results from minor merger with g = 0.1. However in Figured! 
some spin-flip radio sources with apparent jet reorientation angle 15° < 6 < 40° are from 
the minor merger with g = 0.1, this is because of the projection ampliflcation. Comparing 
Figure [Hand Figure |21 one can see that the projection ampliflcation becomes more important 
for larger apparent jet reorientation angle. For mergers with g = 0.3, the maximum intrinsic 
spin flip angle is A < 30° in Figure [21 but in Figure [1] a signiflcant fraction of spin-flip radio 
sources with apparent jet reorientation angle 40° < 6 < 65° are due to these mergers and 
even their contributions to the formation of spin-flip radio sources with 65° < 6 < 90° are 
not negligible because of the projection ampliflcations. The ampliflcation effect becomes 
dominant for the formation of spin-flip radio sources with 65° < 6 < 90°, because Figure [2] 
showed that even for equal-mass mergers g = 1 the probability to have spin flip angle larger 
than 60° is nearly zero. 

It has been suggested that the s pins of SMBHs may be partially aligned v i a interactions 



with a gas disk or via GR precession (iBogdanovic et al.ll2009 



Dotti et al. 2009; Kesden et al. 



2010l ) prior to the binary coalescence. To understand the effects of spin alignment on the 
distribution of spin-flip angle, we ran test simulations assuming that the black hole spins are 
initially partially aligned within 30° relative to the orbital angular momentum of SMBHB 
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(also assumed by iLousto et al.l l2010al ). The corresponding results are shown in the right 
panel of Figure [2l We did not consider further alignment within 10°, because such extreme 
alignment hardly ever produces an apparent spin-flip angle greater than 15°, the minimum 
angle required to identify spin flip in radio sources. When spin alignment is taken into 
account, the peaks of the distribution functions move systematically to lower A and the 
FWHMs shrink significantly. The maximum intrinsic spin- flip angle, ~ 30°, is achieved 
when q = 1. These results imply that in the context of spin-flip model, radio sources with 
spin-flip angles much greater than ~ 30° are strong evidences against spin alignment during 
SMBHB coalescence. 

While the spin axis of jet-ejecting SMBHs reorient during the merging of SMBHBs, the 
post-merger SMBHs get kicked due to the anisotropic radiation of gravitational wave. In the 
upper panel of Figure [31 we give the distribution of the recoiling velocity V in spin- flip radio 
sources for black holes with random initial spin orientation and typical mass ratios q = 0.3, 
0.7, and 1, which are, respectively, the typical mass ratios in spin-flip radio sources with 
apparent jet reorientation angles 15° < 5 < 40°, 40° < 5 < 65°, and 65° < 5 < 90°. Our 
results show that mergers of SMBHBs with larger mass ratios could lead to larger recoiling 
velocities with broader distributions. For major mergers of g = 0.3, the recoiling velocities 
show a distribution ranging from nearly km s~^ to larger than 1200 km s~^ with a most 
possible velocity 600kms~^. For major mergers of q=0.7 or q=l, the recoiling velocities 
V show broader distributions ranging from nearly to larger than 2100 km s~^. The most 
probable kick velocity is about 900 km s~^, while the typical probability for a post-merger 
SMBH to have recoiling velocity larger than 1000 km s~^ is nearly 50 percent. In Figure El 
we also give the results for minor mergers with q = 0.1, which will lead to small spin flip 
of A < 16°. The kicked SMBHs for q = 0.1 have with a very narrow velocity range with 
FWHM(V) ^ 210 km s~^ and the most probable kick velocity is 150 km s~^. Therefore, our 
results suggest that the post-merger SMBHs in spin-flip radio sources with larger apparent 
jet reorientation angle should have larger recoiling velocities, in particular for spin-flip radio 
sources with 15° < 6 < 40° the typical recoiling velocity is about 600 km s~^, while for 
spin-flip radio sources with 65° < 6 < 90° the recoiling velocity is typically about 900 km s~^ 
or larger. If pre-merger spin alignment is considered, the typical spin-flip angle would be 
smaller than 30° according to the right panel of Figure [21 The typical recoiling velocity 
would be significantly smaller, as the thin dashed line shows in the upper panel of Figure [31 

In addition to the magnitude, it is also important to know the direction of the recoiling 
velocity of a kicked SMBH. In the lower panel of Figure [31 we give the distribution of the 
angle a of recoiling velocity relative to the spin axis of post-merger SMBH, or the orientation 
of the active jets. Our results suggest that the kicked SMBHs are inclined to move along 
the active jets and likely to avoid moving in the SMBHB orbital plane. Earlier statistical 
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studies by lLousto et al.l ( 12 01 Obi ) on the spin properties of merger remnants showed the similar 
trend and attributed the ahgnment to that both final spin and recoiling velocity are biased 
toward (counter-)aligning with the orbital angular momentum of the pre-merger SMBHB. 
This explanation implies that if the spin axes of SMBHs could align with the orientation 
of orbital angular momentum prior to coalescence, the bias would be even stronger and the 
correlation between the orientations of final spin and recoiling velocity would be tighter. 
The result from our test simulations with partial spin alignment, which is shown in the lower 
panel of Figure E] as the thin dashed line, supports this implication. For major mergers with 
q > 0.3, most kicked SMBHs move in a direction within 40° with respect to the final spin 
axes or the active jets, which is about or less than the opening angle of the dust torus in 
AGN unification model. Therefore, recoiling SMBHs in most spin-fiip radio sources do not 
move across the dust torus and make direct impact on it, if the active jets are vertical to 
the plane of the dust torus. For equal-mass mergers with q = 1, nearly all kicked SMBHs 
move in the direction a ~ 15° anti-parallel to the final spin angular momenta. These results 
imply that if SMBHBs merge in blazars, the broad emission lines are nearly maximumly 
Doppler-shifted for a given recoiling velocity. For minor merger with q = 0.1, the kicked 
SMBHs move in a broad range of directions and most of them move in the direction of about 
50° relative to the active jets. 



3.4. Line-of-sight velocities and apparent off-center displacements of recoiling 

AGNs 

The projections of spin fiip angles vertical to LOS is the apparent jet reorientation angles 
6 of spin-fiip radio sources and can be given by observations. The projected components V\\ 
of recoiling velocities along LOS can produce Doppler-shifted kinematic offsets of broad 
emission lines with respect to the narrow emission lines of AGNs or to the spectral lines 
of host galaxies, while the other projected component V± perpendicular to LOS leads to 
apparent spatial offsets of AGNs with respect to the center of the host galaxies. Therefore, 
in this section we investigate the Doppler-shifting recoiling velocity V\\ and the apparent 
spatial off-center displacement L± of recoiling SMBHs as a function of the apparent jet 
reorientation angle in spin-flip radio sources. To take into account the deceleration of a 
recoiling SMBH due to dynamical friction against stars in the host galaxy and dark matter 
in the halo, we adopted an 1-D analytic model (see Appendix [A]) to calculate V\\ and L± as 
a function of t^, the elapsed time since a post-merger SMBH gets kicked. 



Armitage fc Natarajaru ( 120021 ) sho wed that t h e mer ger of SMBHB in an AGN would 
truncate the inner accretion disk, and iLiu et al.l ( 120031 ) suggested that the jet formation 
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will be interrupted accordingly and later restart, when the inner disk is refilled after an 
interruption timescale ttr- In this scenario, the total elapsed time after the SMBHB coalesces 
is the sum of ttr and the life time of the active jets t^., i.e. te = + ^tr- In X-shaped radio 
sources, one can estimate the radiation ages of active radio lobes using radiation models 
for relativistic electron plasma in the radio lobes, or calculate the dynamic lifetime tdy of 
jets using the typical length of jets and the typical advancing velocity V^s of hot-spots 
in the active radio lobes (e.g. [Alexander fc Leahyl 119871 ). We can use either of th em as the 



typica l lifetime ta of active jets. To calculate ttr, we adapted the equation (10) in iLiu et al. 



(120031), t 



4 X 10^ q^/^{M/10^ 



Mq) yr, for SMBHBs with circular orbitjj and assumed 



viscous parameter of 0.1 and an opening angle of 0.01 for a standard thin accretion disk. The 
resulting interruption timescale of jet formation for M < 10^ M0 is and much shorter than 
the typical lifetime of the active jets in a giant radio source (GRS, typical jet size ~ 100 kpc 
and age ~ 3 x lO^yr, e.g. §|1]). However, in a compact steep-spectrum (CSS, typical jet size 
~ 10 kpc and age ~ lO^yr) or a gigahertz p eaked-spectrum (GPS, typical jet size < 1 kpc 
and age < lO^yr) radio source ( 10'Dealll998l ). the interruption time Ur become comparabl e 
or greater than the lifetime of the active jets (assuming Vhs — 0.3c, lOwsianik et al.lll999l ). 
thus not negligible in the calculation of tg. 

Radio sources are classified i nto radio galaxies and quasars. In the unification model 
of AGNs (jUrry fc Padovanil Il995l ). radio quasars are those radio galaxies that the angles 
between the relativistic jets and LOS are small so that the accretion disks and the outer 
dust tori are nearly face on. Therefore, in quasars the broad emission line regions can be 
observed directly. However, in radio galaxies, the accretion disks and dust tori are nearly 
edge-on and the broad emission line regions are obscured by dust tori. The division viewing 
angle between radio quasars and galaxies is about 44.4°. In the following, we calculate the 
Doppler-shifting recoiling velocities and apparent spatial offsets for spin-flip radio galaxies 
with 44.4° < W2 < 90° and for spin-flip quasars with 10° < W2 < 44.4°, respectively. We do 
not consider the cases with W2 < 10° because such sources will be identifled as blazars due 
to the close alignment of jets and LOS. 

The upper panels of Figure H] show the cumulative distributions of the LOS velocity V\\ 
in spin-flip quasars. In the calculation, we assumed a typical black hole mass of 10® Mq and 
stellar velocity dispersion of host galaxy 200 km s~^. In the left panel, the solid curve refers to 
the initial kick velocity immediately after SMBHB coalescence, and the other curves refer to 
quasars with increasing jet sizes and ages, namely, GPSs, CSSs, and GRSs. With the increase 



^ In equation (10) of iLiu et all (|2003l ) the factor / = 21.8 and the resonance number n = 5 are adopted 
for binary orbit with eccentricity e = 0.68. While, for a circular binary orbit which we assumped here, f ~ I 
and n — 2. 
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of time since kick and the transition of quasar types from GPSs to CSSs and finally GRSs, the 
deceleration of recoiling SMBH due to the dynamical friction against stars in the host galaxy 
and dark matter in the halo becomes more significant, so both the maximum V\\ and the 
probability for Vji > 10 km s~^ decreases. For GRSs, the bump of the cumulative distribution 
function at V \\ < 70 km s~^ is contribu ted mainly by the SMBHs in the long-term oscillatory 
phase II (IGualandris fc MerrittI l2008l ). In the right panel, we show distribution functions 
of V|| for CSS quasars with three different apparent spin-flip angles, i.e., 15° < 6 < 40°, 
40° < (5 < 65°, and 65° < 5 < 90°. The first distribution function for 15° < 5 < 40° 
resembles the distribution function for all CSSs shown in the left panel, indicating that the 
majority of the spin-flip CSSs have small 6. While the latter two distribution functions 
are significantly higher than the first one, implying that V\\ is generally larger for a spin- 
flip quasar with larger apparent jet reorientation angle. Current spectroscopy with medium 
resolution R ~ 1800, such as that f rom Sloan Digital Sky Survey, can resolve a Doppler- 



shift velocity of Av ^167 km s ^ (lEracleous et a 



20111). A spectral resolution as hig h 



as 46,000, such as HIRES/Keck and UVES/VLT flArav et al.lll998l : iBalashev et al.l 120111 ). 
can even resolve Av ^ 6.5 km s^^. Therefore, with a conservative assumption of a spectral 
resolution of 200 km s~^, the implications from the upper panels of Figure [Hare as follows. 
Among all quasars showing spin-flip signatures, the probability of detecting Doppler-shifted 
broad emission lines drops from as high as 50% for quasars with quenched jet formation and 
activity, to < 20% in GPS and CSS quasars, and to < 2% in GRSs. However, if only quasars 
with 6 > 40° are considered, the detection probability would rise by several times up to an 
order of magnitude. The above results also suggest that spin-flip quasars hardly ever have 
LOS velocities greater than 2000 km s~^. 

In the lower panels of Figure HI we give the cumulative distribution functions of Lj_, the 
apparent off-center displacement, for quasars with different typical jet sizes (left) and with 
different apparent spin-flip angles (right). As the jet size grows and deceleration proceeds, 
the maximum L±_ increase while the probability for L±_ > 1 pc decreases. The probability for 
a spin-flip quasar to have > 2 kpc is negligible, therefore the recoiling quasar is confined 
in its host galaxy. We can also see that for quasars with larger 6, the probability of L± > 1 pc 
is generally greater. To detect an off-center quasar, one needs to observe the angular offset 
between the recoiling SMBH and the center of its host galaxy. Using the images obtained 
with Hubble Space Telescope (HST) of spatial resolution about 100 milli arcsec (mas) and a 
relative astrometric accuracy ~ 12 mas. iBatcheldor et al.l (120101 ) recently detected an angular 
displacement of ~ 100 mas between the SMBH and the photo-center center of host galaxy in 
M87. The detection implies that the relative angular astrometric accuray of a telescope may 
be more important in observing off-center displacements of recoiling SMBHs. The current 
large telescopes, e.g. the Keck II telescope and Very Large Telescope (VLT) equipped with 
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adaptive optics (AO) systems, can achive the highest ang ular resolution ^ 20 mas and relative 



astrometric accuracies ~ 0.2 mas at infrared wavelength (jWizinowich et al. Il2000l : iLu et al 



2010l ). while future astro metric interferomete r missions, such as Gaia, are designed to achieve 



sub-mas resolution (e.g. iPopovic et al.ll201l[ ). At a distance of 100 Mpc (1 Gpc), an angular 
displacement of 0.1 mas corresponds to a physical scale of 0.05 (0.5) pc. The lower panels 
of Figure H] indicates that the probability with large telescopes with AO systems or Gaia 
detecting off-centered quasars among the spin-flip quasar sample is > 50% (> 30 %) for 0.05 
(0.5) pc. In the GPS subsample, especially, the detection probability is 95% both for 0.05 
pc and 0.5 pc. 

Figure |5] shows the same distribution functions as in Figure IH but for spin-flip galaxies. 
The distribution functions show similar dependence on typical jet size and apparent spin-flip 
angle as that in spin-flip quasars. However, due to the projection effect, the maximum and 
mean values of V\\ for spin-flip galaxies are smaller than those for spin-flip quasars, and the 
maximum and mean values of L± are slightly greater. 

The above results are derived assuming a black hole mass of 10® Mq. The corresponding 
kick velocity for a recoiling SMBH to escape to the effective radius of the host galaxy is about 
1600 km s~^. For systems with larger (smaller) black holes, the escape velocity becomes 
higher (lower). As a result, given t^, the probability of detecting either Doppler-shifted 
broad emission lines or off-center AGN displacement would decrease (rise). We did not 
consider spin alignment in this subsection, which (if taken into account) would reduce 
and the possibilities of detecting appreciable V\\ or L^. Since SMBHB coalescence with spin 
alignment hardly ever produces any radio source with A > 40°, radio sources showing large 
apparent spin-flip angles 6 > 40° are consistent with random spin orientation, and for these 
sources our conclusions are unchanged. 



3.5. Truncation of dust torus and dust poor AGNs 



The dust torus around central SMBH is an essential component in the unification model 
of AGNs. The inner edge of dust torus, rtor.inj is most probably determined by the dust 
sublimation radius 



'"torj 



> 0.4 



^bol 



lO^^erg s 



/ISOOKy-' 

V ^sub / 



(10) 



(e.g. 



Nenkova et al.ll2008l ). where Lboi is the AGN bolometric luminosity, and Tsub is the dust 
sublimation temperature. The outer radius of d ust torus, rtnr.n nt, is poorly determined but 



likely not much larger than rto 



20rto 



(e.g. 



Elitzur 1120071 ). Inside the dust torus is the 
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broad emission line region, the size of which can be given with the empirical scaling relation 

"AL,(5100A)^°'' 



rBLR = 1.9 X 10" 



pc, 



where La(5100A) is the AGN luminosity at A = 5100A flKaspi et al.lboosh . The dichotomy 
of Type II and Type I quasars is then caused by whether or not LOS toward the broad line 
regions of AGNs is blocked by the dust tori. 



Komossa fc MerrittI (j2008al ) showed that in a large fraction of recoiling quasars the 
matter in dust tori would become unbound to the recoiling SMBHs and the Type II quasars 
obscured by the dust tori would transform to unobscured Type Is. When a SMBH gets a 
kick velocity V, the matter orbiting the SMBH at radius r > becomes unbound and only 
those with r < tk remains bound to the recoiling SMBHs. The critical radius is 



GM 



0.43 ( - 

v: 



M 



(- 



V 



WMq V 103km s 



pc, 



(12) 



where Mq is the solar mass and M is the mass of kicked SMBH after merger. If V is larger 



than a critical value Vcr which satisfies rk 



no 



all the matter in dust torus becomes 



unbound, therefore the recoiling quasar must be "dust free" , and the initial radio galaxy will 
become an off-center quasar. If the recoiling velocity V is not only larger than Vcv but also 



wl 



4.8 X 10'^ 



(- 



M 



V lO^Mp, 



1/2 



ALa(5100A) 
lO^^erg s-i 



-0.35 



km s 



(13) 



a part of the matter in broad line region may also become unbound and the recoiling quasar 
will become a weak emission line object. However, if the recoiling velocity is moderate and 
^ ^ ^dp ~ Kr('"tor,in/'"tor,out)^^^, the Critical radius rk is larger than rtor.in but smaller than 
'"tor,out- Then part of the dust torus of the recoiling AGN will become unbound, and the 
object will become "dust poor". 

For a typical radio source of M = 2 x IQ^Mq, Tg^b = ISOO/T, /tor = ''"tor,out/^tor,in ~ 20, 
and Lboi = lO^^ergs"^ or ALa(5100A) ~ lO^^erg s~^, the critical recoiling velocity for 
the dust torus to be completely truncated and the object to become "dust free" is Va = 
1466 km s~^. Meanwhile, the critical velocity for the dust torus to be truncated partially 
and the object to become "dust poor" is Vdp ~ 300 km s~^. Our calculations show that the 
probabilities for such a typical radio source to have velocity larger than the critical velocity Vcr 
are approximately 0.9% for 15° < 5 < 40°, 11% for 40° < 5 < 65°, and 10% for 65° < 5 < 90°. 
Therefore, only about 10% of spin-flip radio sources with 5 > 40° will become "dust free" 
quasars or broad emission line radio galaxies if broad line regions are present. However, the 
critical velocity Vdp is at least two times smaller than the most-probable velocity in spin-flip 
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radio sources (c.f. Figures [T] and |3]) . Our calculations indicate that the probabilities for a 
typical spin-flip radio source to become "dust poor" are approximately 72% for 15° < 6 < 40°, 
85% for 40° < 6 < 65°, and 85% for 65° < 6 < 90°. Therefore, most spin-flip radio sources 
should be "dust poor". To become a weak emission line quasar, the recoiling velocity in a 
typical spin-flip radio source should be ^ > Vwi — 6.7 x 10^ km s~^, which is larger than the 
maximum possible recoiling velocity produced during the merger of SMBHB with circular 
orbit. Since the bolo metric luminosity Lboi should not greatly exceed the Eddington limit 
~ 10^^(M/10* Mq) erg s~^, equation (I13p imposes a lower limit on V„\ for a given black hole 
mass, which is ~ 9.5 x 10^ f^f{M/10^ Mq)^-^^ km where /boi = Lboi/[ALA(5100A)] 
is the bolometric correction factor typically in the range 5 < /boi < 20 for AGNs. For less 
massive recoiling black holes, this lower limit would be smaller and more easily achievable. 
Therefore, weak-line AGNs, if exist, preferentially reside in AGNs with small black holes 
(M ~ 10® Mq) and high luminosity. 



4. Recoiling SMBHs in a sample of X-shaped radio objects 

In section [3|, in order to investigate the distributions of Vji and L± for spin-flip radio 
sources, we have to specify the value or make assumptions for the distributions of the un- 
knowns, such as q, ai, a2, afin, ta, W2, etc. However, for a real X-shaped radio object, 
constraints on the distributions of several unknowns can be made by observations. For ex- 
ample, according to section [3l2| the formation of X-shaped radio feature requires mass ratio 
q > 0.1 and most probably q > 0.2. The fact that X-shaped sources are radio loud AGNs 
implies ai > 0.9 and agn > 0.9. In this section, we study the conditional probability distri- 
butions of V|| and L± under more observational constraints for a sample of real X-shaped 
radio objects. 



4.1. The sample of X-shaped radio objects 



Mezcua et al.l (120111 ) studied a sample of 29 X-sh aped radio ob jects, which is drawn 



from the list of 100 X-shaped radio source candidates (ICheungI 120071 ). They estimated for 
their sample objects the dynamic ages of the active radio lobes and the central black hole 
masses. The central black hole masses are computed with the empirical relationships either 
between SMBH mass and the stellar velocity dispersion of host galaxy, o r between the width 
of broad emission line and optical luminosity XL^^^^^ at A = 5100A. iLandt et al.l ( 2010l ) 
studied the optical spectra of another sample of X-shaped radio objects and estimated the 
viewing angles W2 of the active jets, based on the relationship between the viewing angles and 
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the rest frame equivalent widths of the narrow emission hnes O III A5007 and O II A3727. 
We cross-identified the two samples of X-shaped radio objects and obtained 23 objects to 
constitute our sample of X-shaped radio objects for further study. 

We read the v iewing angles W2 for the 23 X-shaped radio sources from the Figure 1 of 



Landt et all (120101 ) in which they divided the objects into four groups according to 0° < 
W2 < 15°, 15° < W2 < 35°, 35° < ^2 < 60°, and 60° < W2 < 90°. We will follow the 
division o f the objects because of the large uncertainties in the estimations of the viewing 



angles. In iLandt et al.l ( l2010l ). a few of our sample objects were classified into the group of 
60° < W2 < 90°, while upper limits, W2 ^ 70°, were presented in their plot. Therefore, we 
assign these objects in a range of 45° < W2 < 70° in our calculations. 

We measured the apparent jet reorientation angles 6, the angle between the orientations 
of active lobes an d wings, of the 23 X-shaped radio objects from the radio ima ges given by 



Lai fc Raol (120071 ) for 3C192 4C+32 25, 4C +48.29, 1059+169, and 3C223.1, bv lWang et al. 



( l2003l ) for 4C+01.30, and by ICheung] (120071 ) for the rest 17 objects. To do the measurement, 
we have to define the orientations of the active lobes and the wings. The orientation of 
active lobes is defined with the line passing through the center of the host galaxy and the 
hot spots inside the active radio lobes. However, it is very difficult to define the orientation 
of wings because they are very wide, diffusive, and without hot spots inside or central bright 
plasma jets along them. To overcome this problem, we define two orientations of the edges 
of a wing, both of which starting from the galactic center but one along the most furthest 
side (contour) and the other along the closest side (contour) of the wing with respect to the 
active jets. With the two orientations, we then measured two angles Si and 62 between the 
orientations of active lobes and the two edges of wings, respectively, and took the mean value 
of 61 and 62 as the value of apparent jet reorientation angle 6. Because of the limited image 
quality, the obtained values for 6 have very large uncertainties. We take 15° as the typical 
uncertainty for the measurement of 6, but for some objects the uncertainties are much larger. 
In particular, the uncertainties of 5 for the objects J0813+4347, J1043+3131, Jllll+4050, 
and J1210-I-1121 are extremely large, therefore these objects will not be included in our 
further calculations. Finally, our sample of X-shaped radio sources consists of 19 objects 
which are listed in Table [U 

In TableHl Column 1 lists the lAU names of the objects based on the J2000.0 coordinates 
and other common catalog names, and Column 2 gives the spectroscopic identification as 
well as redshift, G for radio galaxies and Q for radio quasars. Columns 3, 4, 6, 7, and 8 
are, respectively, the masses of central SMBHs, the optical luminosities of AGNs at 5100A, 
the viewing angles W2 and the dynamic ages of active lobes, and the references (\4s = 
0.1c for GRSs is assumed in the reference literatures). In Column 5, we give the angles 6 
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between the active lobes and wings and the uncertainties measured in this paper. Column 9 
gives the critical velocities in unit of 100 km s~^ for AGNs to become dust-free, calculated 
according to = rtor.in and assuming /boi = 10. Column 10 gives the critical velocities 
for AGNs to become dust-poor, calculated according to = rtor,out and assuming that 
/tor = 20. Column 11 is the viewing angles wi of the wings. Because wi cannot be constrained 
observationally and could be any value in the range of 20° < wi < 90°, we considered three 
possible bins, namely (20°, 40°), (40°, 60°), and (60°, 90°). 



4.2. Spectral kinematic offsets and apparent spatial displacements 

To compute the projected recoiling velocities along and vertical to LOS using Equa- 
tions (EJ-dH]), we have to know the mass ratio q, spin vectors ai and a2, and the LOS vector 
eg for each object. From observations we can estimate only the quantities wi, W2, and 6, 
and may assume ai > 0.9 and Ofin > 0.9 as in section [3J For the mass ratio q, the results 
given in section IX^ suggest that g > 0.1 and most probably q > q^in ~ 0.2. Therefore, 
we may put a lower limit for q in our simulations. Our simulations showed that different 
lower limits of mass ratio q do not change the conclusions significantly. We do not have 
any further knowledge of q, ai, a2 and Cg. Therefore, we assumed possible distributions of 
the quantities and used Monte Carlo simulations to generate their values. For simplicity, 
we assumed flat probability distributions of the parameters q with gmin < Q' < 1, with 
0.9 < ai < 0.998, and a2 with < 02 < 0.998. We also assumed random orientations of spin 
vectors ai, a2, and the LOS vector eg. Given ai, a2, and eg, we calculated wi, W2, and 6 
according to the relevant formula in section |2l and then we used the observed Wi, W2, and 
6 given in Tabled] to select viable solutions for (ai, a2, eg). Finally, given q, M, and ta, 
we derive te according to the scheme explained in § 13. 4[ and then compute the decelerated 
LOS velocity V\\ and apparent off-center displacement L±, using the 1-D analytical model 
described in Appendix [X] 

As in section[31 Monte Carlo simulations are carried out for each of the 19 X-shaped radio 
objects in our sample. We derived the probability distributions of Doppler-shifting recoiling 
velocity and apparent spatial off-center displacements in the objects, based on the 10^-10^ 
sets of trial solutions to the equations and excluding those which are inconsistent with the 
observational constraints. The final results are obtained with about 10^-10^ trial solutions. 
In our simulations for the objects JllOl+1640, J1210-0341, J1327-0203, and J1444+4147, 
no trial solution of the equations ([3])- ([8]) is found to fit all the observational conditions. This 
may imply that either our assumptions for the physical quantities are invalid or the spin-flip 
is not applicable for these objects. For the rest 15 objects in our sample, final solutions to 



- 20 - 



the equations can be obtained at least for one of the three bins of viewing angle, namely 
20° < u;i < 40°, 40° < u;i < 60°, and 60° < Wi < 90°. The cumulative distributions 
of the LOS velocity V\\ and the apparent offset-center displacement for the 15 objects 
are given in Figs. |MH1 Based on the cumulative distribution functions, we then derive the 
probability of displaying Doppler-shifted broad emission lines, P„j|, and off-center spatial 
displacements, Pe,±- These two probabilities are given, respectively, in Columns 12 and 13 
of Table [TJ The two values in Column 12, from left to right, correspond to V\\ > 500 km s~^ 
and Vji > 200 km s~^, while the two values in Column 13 correspond to angular sizes 9± > 10 
mas and > 0.1 mas. 

The probability Pt,j| for our sample of X-shaped radio sources is typically smaller than 
10%, except for J0941-0143, 4C+01.30, J1348+4411, and J1614+2817. While P^,|| for the for- 
mer three object are still smaller than 20%, the last object, J1348+4411, shows an appreciable 
probability 40%) of having V\\ > 200 km s~^. In the case of Pe,^, for 6± > 10 mas, two 
objects have these probabilities greater than 20%, i.e., J0941— 0143 and J1348+4411. When 
> 0.1 mas is considered, all the 15 objects with viable solutions have Pg^± > 20%. In par- 
ticular, four objects, 4C+32.25, J0941-0143, J1348+4411, and J1614+2817, have Pe,± > 50% 
when 6± > 0.1 mas. Therefore, J1348-I-4411 would be intriguing object for follow-up ob- 
servations, in the sense that it has an appreciable probability to show both Doppler-shifted 
broads emission lines and off-center spatial displacement. 



4.3. Truncation of dust torus and dust poorness of the sample objects 

Given Vcr and Vdp for each object, we also calculated the probability P^f of an X-shaped 
radio object to have initial recoiling velocity larger than V^r and the probability P^p to have 
recoiling velocity larger than Vdp. The probability Pdf is also the probability for a sample 
object to have type transition from radio galaxies to quasars and to become "dust free". 

We did not present the probability Pdf in Table [H because it is less than 20% for 
most of the objects except for the quasar 4C+01.30 and radio galaxy J1348-I-4411. The 
probability Pdf for 4C-I-01.30 to have recoiling velocity larger than its critical velocity ~ 
700 (1000) km is about 50% (30%), assuming /boi = 20 (5). In particular, for 40° < 
Wi < 60° (so that torus is "edge-on" in the initial configuration) the resulting probability 
for this object to be dust-free is about 34% (60%). This result implies that the AGN could 
be obscured before binary black hole merger, but the SMBH may have left behind the entire 
dust torus due to recoil and now transformed into a quasar. While for J1348-f4411, its low 
critical velocity ~ 560 (800) km s~^ suggests that Pdf ~ 63% (45%), thus the object 
may transform to a quasar or broad emission line radio galaxy. However, the probability 
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for the object to remain as a radio galaxy is 37% (55%) and not negligible. Therefore, 
non-detection of broad emission lines does not imply that the SMBHB spin-flip model is 
completely excluded for this object. 

The probability Pdp of an X-shaped radio object to have recoiling velocity larger than 
Vdp and to become "dust poor" is given in Column (14) of Tabled! assuming that /boi = 10 
and /tor = 20. The upper and lower indices show the upper and lower limits to Pjp when a 
range of bolometric correction factor 5 < /boi < 20 and torus size factor 10 < /tor < 40 are 
considered. Among the 15 objects with valid solutions, 12 have Pdp > 50%, and six of them 
even have Pdp > 80%. The upper limit of Pdp is greater than 70% for almost all objects, 
and even greater than 95% for 4C-I-01.30 and J1348+4411. The poornes s of dust in AGN 



can be observationally tested w ith the aid of IR satellites such as Spitzer f lJiang et al.l 12010 



Hao et al. 


2010; 


Hal et al. 


2011) 



4C+32.25, J0941-0143, 4C+01.30, J1207+3352, and J1348+4411, would be interesting tar- 
gets for IR follow-up observations. We note that for 3C192, 4C+32.25, and J1207+3352, the 
probabilities for the off-center displacement L± to be smaller than the typical size of dust 
torus are > 90%. The SMBHs in these objects are likely settling back in the galaxy nuclei 
and may become obscured by a new-forming, next-generation dust torus. 



4.4. Comments on individual objects 

Figures [MH] show the cumulative distribution functions of the LOS velocity V\\ and of 
the apparent spatial off-center displacement Lj_ in host galaxies for the 15 objects in our 
sample which have numerical solutions fulfilling observational constraints. Here we give 
detailed comments for the objects one by one. 

JOOOl-0033: a radio galaxy at z=0.247 (Figure [6]). The probability of having Doppler- 
shifted broad emission lines with Vji > 200 km s~^ is < 1.5%, consistent with given 
the uncertainties in our model assumptions. The probability of having AGN off- 
center displacement with > 0.1 mas is about 40%. The critical velocity Vdp ranges 
from 450 to 1300 km s^^, when the uncertainties in /boi and /tor are considered. The 
corresponding probability of being "dust-poor" ranges from 76% to 12%. The critical 
velocity Vcr ranges from 2800 to 4000 km s~^, resulting in a negligible probability of 
being "dust-free". 

J0049+0059: a radio galaxy at z=0.304 (Figure [6]). The probabilities of detecting Doppler- 
shifted broad emission lines with V\\ > 200 km s~^ and AGN off-center displacement 
with 6± > 10 mas are consistent with 0. The probability of having AGN off-center dis- 
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placement with 6'_l > 0.1 mas is about 40 ~ 50%. The critical velocity Vdp ranges from 
560 to 1600 km s~^, when the uncertainties in /boi and /tor are considered. The corre- 
sponding probability of being "dust-poor" ranges from 64% to (1 ~ 2)%. The critical 
velocity V^r ranges from 3600 to 5100 km s~^, resulting in a negligible probability of 
being "dust-free". 

J0808-F2409 (3C192): a radio galaxy at z=0.060 (Figure [6]). Notice that the distribution 
functions of Vj] and 9_[_ vary with Wi by as large as an order of magnitude. The 
probability of having Doppler-shifted broad lines with V\\ > 200 km s~^ is non-zero, 
but < 5%. The probability of having AGN off-center displacement with 6± > 0.1 
mas is about 50%. The critical velocity Vdp ranges from 240 to 690 km s~^, when the 
uncertainties in /boi and /tor are considered. The corresponding probability of being 
"dust-poor" ranges from approximately 56% to 89%. The critical velocity V^j. ranges 
from 1500 to 2200 km s~^, resulting in a negligible probability of being "dust-free". 

J0831+3219 (4C+32.25): a radio galaxy at z=0.051 (FigureE]). For 60° < wi < 90°, none of 
the 10* trial Monte Carlo simulations produces a solution fulfilling all the observational 
constraints. While for other wi, the probabilities of detecting Doppler-shifted broad 
emission lines with V\\ > 200 km s~^ and AGN off-center displacement with 6j_ > 10 
mas are consistent with 0. The probability of having AGN off-center displacement 
with 6± > 0.1 mas is about 50~63%. The critical velocity Vdp ranges from 200 to 
550 km s~^, when the uncertainties in /boi and /tor are considered. The corresponding 
probability of being "dust-poor" ranges from approximately 66% to 87%. The critical 
velocity V^r ranges from 2000 to 5500 km s~^, resulting in a negligible probability of 
being "dust-free". 

J0859-0433: a radio galaxy at z=0.356 (FigureE]). The probability of V\\ > 200 km s"^ is 

< 8%. The probabilities of have AGN off-center displacement with 6*^ > 10 mas is 

< 6%, while the probability for 6± > 0.1 mas ranges from 46% to 60%. We do not have 
data about the bolometric luminosity for this object thus cannot give an estimation to 
the probability of being "dust-poor" or "dust-free" . 

J0924+4233: a radio galaxy at z=0.227 (FigureE]). For 60° < wi < 90°, none of the 10* 
trial Monte Carlo simulations gives a solution fulfilling all the observational constraints. 
While for other wi, the probability of having Doppler-shifted broad emission lines with 
V^l > 200 km s~^ is consistent with 0, given the uncertainties in our model assumptions. 
The probability of having AGN off-center displacement with 9± > 10 mas is non-zero, 
but < 8%, while the probability for 6± > 0.1 mas ranges from 41% to 54%. The 
critical velocity Vdp ranges from 350 to 980 km s~^, when the uncertainties in /boi and 
/tor are considered. The corresponding probability of being "dust-poor" ranges from 
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approximately 30% to 83%. The critical velocity Vcr ranges from 2200 to 3100 km s ^, 
resulting in a negligible probability of being "dust-free" . 

J0941-0143: a radio galaxy at z=0.384 (Figure[7]). The probability of having Doppler-shifted 
broad emission lines with Vji > 200 km s~^ is < 15%. The probabilities of having AGN 
off-center displacement with 6'^ > 10 mas is about 30% and the probability for 6± > 0.1 
mas is about 57%. The critical velocity Vjp ranges from 190 to 550 km s~^, when the 
uncertainties in /boi and /tor are considered. The corresponding probability of being 
"dust-poor" ranges from approximately 67% to 91%. The critical velocity Vcr ranges 
from 1200 to 1700 km s~^, and the resulting probability of being "dust-free" is ~ 17%. 

J0941+3944 (3C223.1): a radio galaxy at z=0.107 (Figure [ZD. For 60° < Wi < 90°, none 
of the 10® trial Monte Carlo simulations gives a solution fulfilling all the observational 
constraints. While for other wi, the probabilities of having Doppler-shifted broad 
lines with V\\ > 200 km s~^ and AGN off-center displacement with 6± > 10 mas 
are both consistent with < 8%, given the uncertainties in our model assumptions. For 
9± > 0.1 mas, the probability is about 24%. The critical velocity Vdp ranges from 200 to 
580 km s~^, when the uncertainties in /boi and /tor are considered. The corresponding 
probability of being "dust-poor" ranges from 54% to 89%. The critical velocity Vcr 
ranges from 2000 to 5800 km s~^. The resulting probability of being "dust-free" is 
(0.01 ~ 4.4)% for 20° < wi < 40° and (0.08 ~ 12)% for 40° < wi < 60°. 

J1005+1154: a radio galaxy at z=0.166 (Figure [7]). For 60° < Wi < 90°, none of the 
10® trial Monte Carlo simulations produces a solution satisfying all the observational 
constraints. While for other wi, the probabilities of showing Doppler-shifted broad 
emission lines with V\\ > 200 km s~^ and AGN off-center displacement with 6± > 10 
mas are consistent with 0, given the uncertainties in our model assumptions. For 
6± > 0.1 mas, the probability ranges from 23% to 42%. The critical velocity Vdp ranges 
from 350 to 1000 km s~^, when the uncertainties in /boi and /tor are considered. The 
corresponding probability of being "dust-poor" ranges from 25% to 76%. The critical 
velocity V^r ranges from 2200 to 3200 km s~^, resulting in a negligible probability of 
being "dust-free". 

J1020+4831 (4C+48.29): a radio galaxy at z=0.053 (Figure!?]). The probability of detecting 
Doppler-shifted broad emission lines with Vji > 200 km s~^ is consistent with 0, given 
the uncertainties in our model assumptions. The probability of detecting AGN off- 
center displacement with 6± > 10 mas is non-zero, but < 7%. For 6*^ > 0.1 mas, the 
probability is about 45%. The critical velocity Vdp ranges from 220 to 620 km s~^, 
when the uncertainties in /boi and /tor are considered. The corresponding probability 
of being "dust-poor" ranges from 58% to 87%. The critical velocity V^r ranges from 
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1400 to 2200 km s ^. The resulting probability of being "dust-free" is (0 ~ 11)% for 
20° < u;i < 40°, (0 ~ 7.9)% for 40° < Wi < 60°, and (0 ~ 8.2)% for 60° < Wi < 90°. 

J1130+0058 (4C+01.30): a radio quasar at z=0.132 (Figure E]). For 60° < u;i < 90°, 
none of the 10^ trial Monte Carlo simulations produces a solution satisfying all the 
observational constraints. While for other wi, the probability of Vji > 500 km s"^ 
is < 11%, but the probability of Vji > 200 km s~^ ranges from 15% to 20%. The 
probability of detecting AGN off-center displacement is non-zero, ranging from (16 ~ 
24)% in the case 20° < wi < 40° to (25 ~ 33)% in the case 40° < wi < 60°. The 
critical velocity Vdp ranges from 700 to 990 km s~^, when the uncertainties in /boi and 
/tor are considered. The corresponding probability of being "dust-poor" ranges from 
80% to 96%. The critical velocity Vcr ranges from 1100 to 3100 km s~^. The resulting 
probability of being "dust-free" is (25 ~ 50)% for 20° < < 40° and (34 ~ 59)% for 



40° < Wi < 60°. It is interesting to note that IWang et al.l fl2003[ ) reported a possible 
large offset, 3.5 ", in this object between the radio core and the photometric center 
of the host galaxy. Our Monte Carlo simulations suggest that the observation would 
require a recoiling velocity of V > 4000 km s~^ at a confidence level of 99.9%. Such 
high recoiling velocity is extremely difficult to achieve by asymmetric gravitational 
wave radiation of circular S MBHB orbit. This implies that either the two SMBHs 



merge along hyperbolic orbit (IHealy et al.ll2009l ) or a mechanism other than recoiling 



SMBH is needed to produce the large displacement. 

J1140-I-1057: a radio galaxy at z= 0.081 (Figure |8]). The probability of having Doppler- 
shifted broad emission lines with V\\ > 200 km s~^ is generally < 8%. The probability 
of having AGN off-center displacement with 6± > 10 mas is about 18%, increasing to 
about 50% for 6*^ > 0.1 mas. The critical velocity Vdp ranges from 250 to 700 km s~^, 
when the uncertainties in /boi and /tor are considered. The corresponding probability 
of being "dust-poor" ranges from 55% to 87%. The critical velocity V^r ranges from 
1600 to 2200 km s~^, resulting in a negligible probability of being "dust-free". 

J1207+3352: a radio galaxy at z=0.079 (Figure [8]). For 60° < Wi < 90°, none of the 
10^ trial Monte Carlo simulations produces a solution satisfying all the observational 
constraints. While for other wi, the probabilities of showing Doppler-shifted broad 
emission lines with V\\ > 200 km s~^ and AGN off-center displacement with 6± > 
10 mas are both < 15%. For 6± > 0.1 mas, the probability is about 43%. The 
critical velocity Vdp ranges from 160 to 440 km s~^, when the uncertainties in /boi and 
/tor are considered. The corresponding probability of being "dust-poor" ranges from 
approximately 67% to 92%, The critical velocity V^r ranges from 980 to 1400 km s~^. 
The resulting probability of being "dust-free" is (2.7 ~ 17)% for 20° < wi < 40° and 
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(8.6 ~ 32)% for 40° < u;i < 60°. 

J1348+4411: a radio galaxy at z=0.267 (Figure [8]). The probability of V\\ > 500 km 
is about 13%, but the probability of Vji > 200 km increases steeply to about 35%. 
The probability of having AGN off-center displacement is ranging from about 52% in 
the case of > 10 mas to about 63% in the case of 6± > 0.1 mas. The critical 
velocity Vdp ranges from 89 to 250 km s~^, when the uncertainties in /boi and /tor are 
considered. The corresponding probability of being "dust-poor" ranges from 85% to 
97%. The critical velocity Vcr ranges from 560 to 800 km s~^. For all wi, the probability 
of being "dust-free" and become type I AGN is (45 ~ 64)%. The lack of broad line in 
this object implies that V < Kr- If ^ < Kr is required, our Monte Carlo simulations 
suggests that -P(Vj| > 200 km s"^) and P{9± > 0.1 mas) both decrease to < 10%, but 
Pdp is still in the range (70 ~ 90)%. Alternatively, it is possible that this ra dio galaxy 



intrin sically lacks broad emission lines like in some AGNs and quasars (e.g. iFan et al. 



1999|). 



J1614-I-2817: a radio galaxy at z=0.108 (Figure |8]). The probabilities of having Doppler- 
shifted broad emission lines with Vji > 200 km s~^ and AGN off-center displacement 
with 6*^ > 10 mas are non-zero, but < 20%. For 6± > 0.1 mas, the probability 
increases to about 65%. The critical velocity Vdp ranges from 610 to 1700 km s~^, 
when the uncertainties in /boi and /tor are considered. The corresponding probability 
of being "dust-poor" ranges from nearly 0% to 65%. The critical velocity ranges 
from 3800 to 5400 km s~^, resulting in a negligible probability of being "dust-free". 



5. Discussions and conclusions 

A popular model in the literature suggested that the X-shaped radio sources form be- 
cause of jet reorientation following the spin flip of SMBHs at SMBHB mergers. Numerical 
relativity suggested that the black hole merger not only leads to a spin flip but also to a 
gravitational kick of the post-merger black holes. In this paper, we investigated with Monte 
Carlo simulations the distribution function of the spin-flip angles and recoiling velocities of 
the post-merger SMBHs, and the detectability of recoiling SMBHs, in the forms of Doppler- 
shifted spectral lines and off-center AGNs, in radio sources showing spin-flip signatures of 
jet reorientation. We also collected a sample of 19 X-shaped radio objects and calculated 
for all the sample objects the spectral kinematic offsets Vj], the apparent spatial off-center 
displacements L±, and the degree of dust-torus truncation because of the gravitational kick. 

Because the spin flip and recoiling velocity strongly depend on the parameters q, ai, a2. 
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and eg of SMBHB system, we first investigated tfie distribution of mass ratio q in spin-flip 
radio sources witli given apparent jet reorientation angles, assuming that the spin parameter 
a of jet-ejecting SMBHs has a > 0.9. Our results indicate that to cause significant apparent 
jet reorientation the mass ratio cannot be smaller than a minimum value ~ 0.1. The 
larger the apparent jet reorientation angle is, the larger the minimum mass ratio will be. For 
spin-flip radio sources with apparent jet reorientation angle 15° < S < 40°, gmin ~ 0.05 and 
the most probable mass ratio is g ~ 0.3, while for spin-flip radio sources with 65° < 6 < 90°, 
g^in ^ 0.15 and the most probable mass ratio is g > 0.6. Our results suggest that if 
an observed jet reorientation is indeed due to the spin flip following a SMBHB merger, the 
merger must be major merger with q > 0.2 — 0.3. If the spin-flip model for radio sources of jet 
reorientation could be confirmed, our results suggest that we can measure the distribution of 
SMBHB mass ratio in luminous radio sources and give statistic constraints on the hierarchical 
galaxy formation model and the growth history of SMBHs. 

Having the mass ratio distribution, we investigated the distributions of spin flip angle 
and recoiling velocity in spin-flip radio sources. Our results showed that the larger the 
SMBHB mass ratio is, the larger the the most possible spin flip angle and the recoiling 
velocity will be. For typical major merger mass ratio q = 0.3, q = 0.7 and q = 1.0, the 
most possible spin flip angles are moderate and increase with mass ratio from A = 19° for 
q = 0.3, and 25° for q = 0.7 to 30° for q = 1. The spin flip angle is hardly ever larger than 
A ~ 55°. Our results imply that any apparent jet reorientation angle 6 > 60° is due to 
projection effect. For q = 0.3, the recoiling velocity ranges from about zero to larger than 
1200 km s~^ with a most possible velocity V ~ 600 kms~^. While for q > 0.7, the distribution 
of recoiling velocity is similar and ranges from about to larger than 2100 km s~^ with most 
probable velocity V ~ 900 km s~^. For major mergers with q > 0.3, almost all recoiling 
SMBHs move nearly along the final jet orientation within a angle of 40°. Alignment of spin 
axes with orbital angular momentum of SMBHB prior to the black hole coalescence would 
significantly reduce the recoiling velocity and spin-flip angle, but tighten the correlation 
between the orientations of final spin and recoiling velocity. Our results suggest that to 
efficiently detect recoiling SMHBs with Doppler shifted spectral broad emission lines, one 
should preferentially observe nearly face-on AGNs, e.g. radio loud quasars or blazars among 
which the broad emission lines are nearly maximumly Doppler-shifted for a given recoiling 
velocity. This conclusion is also consistent with present observations in the literature that 
all candidates for recoiling SMBHs are bright q uasars and the measu red Doppler-shifting 



recoiling velocity is preferentially very larger (e.g. iKomossa et al.ll2008l ). However, to detect 
recoiling SMBHs in radio galaxies, one should try to measure the apparent spatial off-center 
displacements of AGNs relative to the centers of the host galaxies because the recoiling 
SMBHs are moving nearly along the jet axes, which are at large angles with respect to LOS. 
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Even in some quasars, recoiling SMBHs may be discovered by measuring the apparent spatial 
off-center displacements because some radio galaxies and T ype II quasars may transit to Type 



I quasars due to the recoiling oscillations as suggested by iKomossa fc MerrittI ( l2008al ) . For 
minor mergers with q ~ 0.1, both the spin-flip angle and the recoiling velocity narrowly 
distribute with peaks at the most possible angle A ~ 12° and velocity V ~ 150 km s~^. 

In order to use the observed jet reorientation angles in spin-flip radio sources to constrain 
the observable spectral kinematic offsets V\\ and apparent spatial off-center displacements L^, 
we projected the spin-flip angles and the recoiling velocities along and vertical to LOS and 
investigated their correlations. An analytical 1-D model is used to calculate the deceleration 
of a recoiling SMBH due to dynamical friction, and the elapsed time since the black hole 
gets kicked is estimated from the size of the radio morphology. In general, Vji in spin-flip 
radio sources is smaller than 2000 km s~^ and Lx_ hardly ever exceeds 2 kpc. However, we 
should note that the conclusion is valid only for quasars with SMBHBs of circular orbits. 
For two SMBHs merging along hyper bolic orbits, it is p ossible to have a recoiling velocity 



significantly larger than 2000 km s ^ (IHealy et al.ll2009l ). Our results indicate that in radio 



sources showing spin-flip signatures of jet reorientation, detectable Doppler-shifted broad 
emission lines are mainly contributed by GPS and CSS sources with large apparent spin-flip 
angles. The detection rate of Doppler-shifted broad emission lines in a sample of GPSs/CSSs 
with 5 > 40° could be as high 50%. For GRSs, deceleration due to dynamical friction has 
significantly reduced V\\ , so the probability of detecting Doppler-shifted broad emission lines is 
typically < 10%. Detecting AGN off-center displacement requires astrometric observations 
with high angular resolution. For current telescopes with limit resolution up to 10 mas, 
detection of off-center displacement would be limited to the GPS/CSS sources within a 
distance of 100 Mpc. However, future interferometer missions with sub-mas resolution or 
higher, such as Gaia, would allow detection of off-center AGNs in GPS, CSS, and GRSs at 
a much larger distance. Our results suggest that to efficiently detect the recoiling SMBHs 
in spin-flip radio sources, one should select a sample of GPS/CSS sources with apparent jet 
reorientation angles 6 > 40°. 

When the SMBH in an AGNs gets kicked, the outer part of the dust torus in AGN 
unification model may become unbound and the dust torus is truncated, leading to the 
transition of the AGN from dust rich to "dust poor" or even to "dust free". When a 
FRII radio galaxy or an obscured quasar becomes "dust free", it should transform to an 
unobscured broad emission line quasar. We calculated the minimum recoiling velocities for 
partial or full truncation of dust torus and estimated their detectability in typical spin- 
flip radio sources with different apparent jet reorientation angles 6. Our results suggest that 
about 11% spin-flip radio galaxies with 6 > 40° and about 1% with 6 < 40° will become "dust 
free" and broad emission line AGNs, if broad line regions are present. Our results imply that 
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quasars with larger bolometric luminosity and smaller black hole mass should have higher 
probability to become "dust free". The velocity for dust torus to be partially truncated 
is typically lower than the most-probable recoiling velocity, therefore most X-shaped radio 
sources should be "dust poor" . Quasars may become weak emission line objects, if the broad 
emission regions are partially truncated due to gravitational recoiling. This is possible for 
a quasar with small black hole mass (~ 10^ Mq) radiating at the Eddington limit. We 
note that without gas replenishing the accretion disks, the dust-free and weak-line quasars 
formed by recoiling SMBHs are relatively short-living (< 1 Myr), because the combination 
of high recoiling velocity (> 10'^ km s"^), small black hole mass (< 10^ Mq), and large 
bolometric luminosity (> 10^^ erg s~^) will significantl y reduce the lifetimes of the accretion 



disks around the recoiling black holes (e.g. iLoebl 120071 ). This argument implies that among 



X-shaped radio sources whose host galaxies are typically gas-poor elliptical, the ones with 
long active jets (> 100 kpc or > 3 x 10^ yr for a typical advancing speed 0.1c of hot spot in 
FRII radio galaxies) and small black hole masses (< 10^ M©) are unlikely to host dust-free 
and weak-line quasars. 

We then selected a sample of 19 X-shaped radio objects, consisting of 18 X-shaped radio 
galaxies and 1 quasar, and for each sample object we calculated the distribution functions 
of the LOS recoiling velocity V\\ and the apparent spatial off-center displacement L±, as 
well as the critical recoiling velocities Vcr and Vdp for partial and complete truncation of 
dust torus. To do the calculations, we measured the apparent jet reorientation angles from 
the radio images and collected the central SMBH masses, the optical luminosities of AGNs, 
the viewing angles and dynamic ages of the active lobes from the literature. Our results 
showed that most X-shaped radio objects in the sample have negligible Vji {Pv,\\ < 10%), 
except for J0941-0143, 4C+01.30, J1348+4411, and J1614+2817. All the 15 sample objects 
with valid Monte Carlo solutions have significant probabilities (> 20%) to show AGN off- 
center displacements when a spatial resolution of 0.1 mas is considered. In particular, the 
radio galaxy J1348+4411 have appreciable probability to show both Doppler-shifted broad 
emission lines {Pv,\\ — 40%) and off-center spatial displacement (-Pe.x — 55%), therefore 
would be intriguing target for follow-up observations. Our results also showed that the 
sample X-shaped radio objects have high probabilities to be "dust-poor" . The probabilities 
for X-shaped radio galaxy J1348-I-4411 and quasar 4C+01.30 are always higher than 80%, 
despite the uncertainties in torus size and bolometric correction factor. Because of the 
high probability of being both dust-poor and offset from galaxy center, these two objects 
would be interesting targets for IR follow-up observations to test the dust-poorness. Besides, 
J1348+4411 and 4C-I-01.30 also have the highest probability in our sample to be "dust-free" 
and non-obscured. The fact that 4C-I-01.30 is already a quasar implies that this object could 
be an obscured AGN before binary black hole merger, but may have transformed into a 
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quasar after recoil takes place. Detecting the predicted off-center displacement and Doppler 
shifted broad emission lines would strongly support that 4C+01.30 is a recoiling SMBH. 
For J1348+4411, the lack of broad line in this object suggests that the recoiling velocity 
is toward the lower end of the predicted distribution function, but the probability of being 
dust-poor is still higher than 70%. 

The detection of the predicted observational signatures of recoiling SMBHs in X-shaped 
radio sources would give final confirmation of the spin-flip model for X-shaped radio sources. 
While, if none of the predicted signatures is detected in X-shaped radio sources, the SMBHB 
spin-flip model for the formation of X-shaped radio feature would be most likely excluded, 
unless our assumptions for the initial conditions of SMBHB mergers are severely unrealistic. 
However, this does not imply that X-shaped radio sources do not harbor SMBHBs at center, 
because je t reorientation may form due to the interaction of SMBHB and accretion disk 



( iLiu 1120041 ). According to the disk-binary interaction model, the final recoiling velocity of a 
SMBH is small due to the alignment between the black hole spin axes and the orbital angular 
momentum before the merger, therefore neither the spectral kinematic offset nor the spatial 
AGN off-center displacement are as large as those in the spin flip model. 
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Fig. 2. — Distribution functions of intrinsic spin flip angle A for different SMBHB mass 
ratios, assuming random initial spin orientation (left panel) or initial spin alignment within 
30° relative to the binary orbital angular momentum (right panel). Spin- flip angles in the 
mergers with larger mass ratio generally have broader distributions and the equal-mass 
mergers give the largest possible spin flip angle. When the initial spin orientations are 
random (aligned), typical spin flip angles are A ~ 30° (~ 19°) for g = 1, A ~ 25° (~ 17°) 
for q = 0.7, and A ~ 19° (~ 11°) for q = 0.3. While for minor mergers of g = 0.1, the spin 
flip angles narrowly distribute with FWHM(A) ~ 7° (~ 3°) and a peak at A = 12° (~ 4°). 
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Fig. 3. — (Upper) Distribution functions of recoiling velocity V for different black hole 
mass ratios. For major mergers with q > 0.7, the distributions of recoiling velocity are 
almost identical, ranging from to about 2100 km s~^ with most possible recoiling velocity 

V ~ 900 km s~^. For smaller mass ratios, V distributes in narrower ranges and the fraction 
with large recoiling velocity decreases rapidly. Most mergers have recoiling velocity V ~ 
600 km s~^ for q=0.3. for minor mergers q ~ 0.1, recoiling velocity is narrowly peaked at 

V ~ 150 km s~^. (Lower) Distribution functions of the direction angle a of recoiling velocity 
relative to the spin axis of post-merger black hole for different mass ratios q. For major 
mergers of g > 0.3, recoiling SMBHs move preferentially along the spin axis of post-merger 
SMBHs within an angle of 40°. For major mergers q > 0.7, most recoiling SMBHs moves in 
a direction of a > 165°. The thin dashed lines in both panels show the results when partial 
spin alignment is considered. 
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Fig. 4. — Cumulative distribution functions of LOS velocity ( Upper) and apparent off-center 
displacement (Lower), for spin-flip quasars at initial kick (with quenched jet formation) and 
with increasing active jet sizes and ages from GPSs to CSSs and finally GRSs (left) and for 
CSS quasars with different apparent spin-flip angles (right). The fiducial black hole mass is 
10^ M0 and the stellar velocity dispersion of host galaxy is 200 km s~^. 
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Fig. 5.— 



The same as FigurelUbut for spin-flip radio galaxies, which have 44.4° < W2 < 90°. 
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Fig. 6. — Cumulative distribution of LOS velocity Vji (Left) and apparent spatial off- 
center displacements L_i_ vertical to LOS (right) for X-shaped radio objects JOOOl-0033, 
J0049+0059, 3C192, J0831+3219, and J0859-0433. The dots and circles in the left panels 
mark, respectively, the locations of V\\ = 200 km s~^ and 500 km s~^, while those in the right 
panels mark, respectively, the locations of 6*^ = 10 mas and 0.1 mas. 
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Fig. 7.— The same as Figure E but for objects J0924+4333, J0941-0143, 3C223.1, 
J1005+1154, 4C+48.29. 
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Fig. 8. — The same as Figure [6] but for 
J1348+4411, and J1614+2817. 
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objects 4C+01.30, J1140+1057, J1207+3352, 
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A. Deceleration of recoiling SMBH due to dynamical friction 



We set up an 1-D analytic model to calculate the deceleration of a recoiling SMBH 
due to dynamical frictio n against stars and dark matter (DM). We adopted a core-Sersic law 
( Terzic fc Gralianj|2005[) for the stellar density distribution and an NFW profile for DM halo 
( iVolonteri et al.ll2003l ). The normalization of the density is determined by the conditions that 
(1) the mass ra tio between SMBH and host galaxy is M,/Mg = 0.002 and (2) the mass deficit 
(lMerrittll2006l) inside cor e radiu s equals the mass of SMBH. Interested readers are referred 



to IGualandris &: MerrittI (120081 ) for more detailed discussions on the density distribution of 
host galaxy. 

For a recoiling spuermassive black hole with mass M, at galactic centers r = obtaining 
an in itial kick velocity v = Vk, the N-body numerical simulations ( iMerrittl l2006l : iLi et al. 



20 111 ) suggests that the dynamical evolution of recoiling black hole within the host galaxy 
and dark halo can be divided into three phases and is well described by the equations 
incorporating the Chandrasekhar formula for dynamical friction 

dr 

— = V 
dt 



dv 
IE 



GMg{r) 



An\YiKG'^M,p{r) 



erf(X) 



2X 



-X2 



(Al) 
(A2) 



where Mg{r) is the total mass inside the sphere of radius r about the galaxy center, In A is 
the Coulomb logarithm, p(r) is the total density of stars and dark matter, X = v/ \/2a{r). 



and a {r) is the stellar velocity dispersion at r. The N-body simulation given by iMerritt 



(120061 ) shows that when the total mass within the orbit of the oscillating black hole is larger 
than the black hole mass M,, the dynamic evolution of recoiling black hole is at Phase I and 
the fiducial Coulomb logarithm is In A = 2.5. However, when the stellar mass within the 
orbit of black hole become smaller than M,, the dynamic evolution of recoiling black hole 
goes into an long-term oscillatory Phase II during which the Coulor nb logarithm is in the 
range 0.1 < InA < 0.3 with a typical value InA ~ 0.2 (Figure 10 of IGualandris fc Merritt 



20081 ). At a much later time of Phase HI, the black ho le reaches a thermal equilibr ium with 
stars and has a Brownian motion at the galaxy center ( IGualandris fc Merrittll2008l ). For the 
cases which we are interested in, the black hole never goes into Phase HI and we will not 
considered it here. 

Figure |9] shows the resulting evolution of r as a function of time t, where r is normalized 
by the effective radius of the host galaxy, i?e, and t is normalized by the dynamical timescale 
Te at Re- Different trajectories are for different initial kick velocity f^, which is normalized by 
the escape velocity v^sc o f the host galaxy and l i sted in the legend. Our model is essentially 
the same as Model A2 in IGualandris fc MerrittI (120081 ) . except that the ratio M,/Mg in our 
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''there are no solutions in our Monte Carlo calculations with 10* random number experiments 
''This source belongs to X-shaped quasar 
'^no SDSS photometry available 



References. — (ll lCheund bOOTi '). f2) lLal fc RacI bOO?!') . (3') lLandt etai] bOldl ). (4') lMezcua et al.l ll201ll '). (5) IWang et~all bOOSl) 



model is tw o times greater, which result s in a slightly longer dynamical friction timescale 
than that in ICualandris fc MerrittI koO^ . 



Notice that in our simple analytic model, the normafeec? trajectories of r / R(.{t/T(.) and 
v/vcscit/Tf,) are uniquely determined by Vk/vcsc and do not depend on M,. This characteris- 
tics allows us to derive the trajectories of r(t) and v{t) for any kick velocity Vk and black hole 
mass M,, by interpolating among a limited number of normalized trajectories pre-calculated 
for different Vk/v^sc values. The normalized trajectories that we used in our Monte Carlo 
simulations are for Vk/v^sc = 0.01, 0.05, 0.1, 0.2, 0.3, 2.0. 
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Fig. 9. — Trajectories of recoiling SMBH with different initial kick ve locity. Because of the 
logarithmic scale, the oscillatory phase II ( iGualandris fc Merrittll2008l ) is clearly resolved in 
the plot when Vk/vesc = 0.3 and 0.7. 



